Elastic properties and phase stability of AgBr under pressure. 
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Ab initio calculations have been used to derive the elastic constants, equation of state and free 
energy of NaCl, KOH and CsCl-type structures of AgBr. The elastic constants have been derived 
by the stress-strain relation. The shear elastic constants C44 has proved to be exceptionally soft. 
We have found that at temperature T=0 the phase transition from the low pressure NaCl-type to 
the high-pressure CsCl-type structure goes over a monoclinic KOH-type structure. The free energy 
comparison indicates that the monoclinic phase ranges from 8 to 35 GPa. 

PACS numbers: 05.70.F, 64.60, 64.70, 71.15.Mb 



Silver halide crystals are of considerable interest, since 
they exhibit deviation from the ideal ionic character due 
to the presence of the d-electrons in Ag + ions. Numer- 
ous studies of these crystals have been undertaken in 
the past years using various forms of interaction poten- 
tials including three-body interactions and van der Waals 
forces)^* 2 - as well as the ab initio total energy pseudopo- 
tential calculations^ There are also many experimental 
works studying various aspects of these solids 4 ri±2r22 
The results of studies using phenomenological potentials 
generally show good agreement with experiment^ This 
studies use sophisticated schemes involving three body 
interactions and van der Waals potentials. The ab initio 
results for transition pressure seems to be less accurate^. 
There is one ab initio study of phase stability of silver 
bromide AgBr.™ Nevertheless, the transition mechanism 
is still unclear for this compound. Several models have 
been proposed^— ii2r— but the transition mechanism re- 
mains still an open problem. The ab initio study by 
Nunes, Allen and Martins ji» which investigated the en- 
ergy of several candidate structures, proposed an inter- 
mediate trigonal phase with space group P3i21. In con- 
trary recent experimental work of Hull and Keeni^ sug- 
gested an existence of an intermediate phase with a mon- 
oclinic, KOH-type structure and shear instability as the 
transition mechanism. The silver bromide AgBr is sta- 
ble at ambient temperature and pressure in the rocksalt 
structure Fm3m. In ambient temperature and high pres- 
sures it has a second stable phase with CsCl structure 
Pm3m. The phase transition in between these phases 
has been reported^— to occur at P f =8.3 GPa. In this 
paper we are studying elastic properties and phase sta- 
bility of the AgBr as a function of pressure. For NaCl and 
CsCl structures we have scanned pressures in the interval 
P=0-200 GPa, and calculated the free energies and elas- 
tic constants. In the range of P=0-40 GPa we have also 
studied the phase stability of the free, non-constrained 
by any symmetry, AgBr crystal , searching for any inter- 
mediate phases as well as general transition mechanism. 
The non-constrained crystal was treated as having PI 
symmetry, however, its unit cell content was not allowed 
to change. Hence, phase transitions to phases with mul- 
tiplication unit cells were not taken into account. 



I. METHODOLOGY 

The calculations have been performed using Density 
Functional Theory (DFT) approach with ultrasoft pseu- 
dopotentials and Generalised Gradient Approximation 
(GGA)^-^ We have used CASTEPSi implementation 
of this method and pseudopotentials provided with this 
package. The elastic constants of the NaCl phase have 
been also verified by the independent calculation using 
VASP 2 - 2 -. The electronic minimisation used density mix- 
ing scheme and pseudopotentials parametrised in the re- 
ciprocal space. The summation over the Brillouin zone 
has been performed with weighted summation over wave 
vectors generated by Monkhorst-Pack schemed All com- 
putations were done in lxlxl supercell with 8 or 2 or 
8 atoms for NaCl or CsCl or KOH structures, respec- 
tively. Some tests for NaCl has been carried out on the 
primitive, others on the face-centered cubic unit cell. 

Due to high demands of the accuracy of the free en- 
ergy, we have conducted extensive tests of convergence 
of the energy and pressure with respect to energy cutoff 
and Brillouin zone grid spacing parameters. Thus, we 
have performed a series of calculations with cutoff ener- 
gies in the range of E cut =300 - 450eV and k-space grid 
spacing in the range Ak=0.02 - 0.1 A -1 . In the light 
of obtained results we have decided to use in our com- 
putations a 0.4 A -1 grid spacing and 340eV energy cut- 
off. The chosen grid corresponds to 32 wave-vectors for 
NaCl structure and 256 wave- vectors for CsCl structure. 
Such a large number of wave-vectors made calculations 
rather time-consuming. But since we are comparing en- 
ergies of different structures with different unit cell vol- 
umes and Brillouin zone shapes, sizes and k-sampling, we 
need a very good convergence of energy. At zero pressure 
the lattice constants of the optimised NaCl structure is 
a= 5.7859 A, while the experimental lattice constants at 
P=0 GPa and T=300 K is A=5.7745 A. 



II. FREE ENERGY AND EQUATION OF STATE 

To investigate the phase stability of the crystal and 
to estimate the phase transition pressure, one must eval- 
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FIG. 1. Lower part: Free energy of AgBr in NaCl (+, solid 
line) and CsCl (x, dashed line) structures. Upper part: 
Free energy difference FcsCi — Gn^c; between CsCl and NaCl 
structures as a function of pressure. The vertical arrow de- 
notes predicted phase transition pressure. 
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FIG. 2. Variation of relative volumes with pressure for 
NaCl(+) and CsCl (x) structure. Experimental points (•, A) 
are taken from Bridgman 5 >>£ and Vaidya^ respectively. The 
dashed line shows results of theoretical study of Gupta and 
Singh-'. The lines are present calculations. 



uate the free energy of all involved phases and choose 
the one with the lowest energy as stable. Since in ab 
initio calculations the temperature T=0, the free energy 
is: F = E + PV. We have performed calculations of 
the crystal free energy F for several pressures from 
to 50 GPa for both NaCl and CsCl structures. Fig [1] 
shows the results of these computations. In the lower 
part the absolute values of free energies per atom are 
plotted together with cubic polynomials fitted to these 
data. In the upper part the difference of FcsCi — Fpj a ci 
is plotted against pressure. The difference changes sign 
at the pressure Pt=27 GPa (arrow). We see that our 
calculations correctly predict that the rocksalt structure 
is more stable at low pressure and that at higher pres- 
sures the stability passes to CsCl structure . The large 
difference between experimental Pt=8.3 GPa and calcu- 
lated P t =27 GPa transition pressure could be attributed 
to a few factors. The main reason could be related to 
the appearance of the intermediate KOH-type structure, 
discussed below. Here, we mention other possible fac- 
tors. First, the experimental data 5 -'- are measured only 
at ambient temperature, while the calculations were per- 
formed at T=0. At finite temperature the entropy con- 
tribution to the Gibbs free energy G = E + PV — TS can 
also be different for the two structures, since the clos- 
est coordination shells contain different number of ions. 
Another reason could be related to the fact that the van 
der Waals interaction is neglected in DFT calculations, 
and that leads to under-binding of the crystal. This ef- 
fect has already been noticed for silver iodide Agl crystal, 
where DFT predicts the phase transition pressure to be 
10 GPa^ while experimentally it is 4 GPa. 

The data used to calculate free energy provide also the 
equation of state of the crystals. Fig. [5] shows the calcu- 
lated equation of state for NaCl and CsCl structures as 
well as the experimental data measured by Bridgman and 



Vaidya^iiS We have also included in figure [2] the results 
of phenomenological study of Gupta and Singhi^ This 
work was able to predict the transition pressure much 
better, but it fails to reproduce the volume compression 
of the phase transition. Present calculations show excel- 
lent agreement with experimental data for the low pres- 
sure range up to 6 GPa. In this interval the relative error 
in V(P)/V(0) ratio is below 0.5%. It should be noted that 
the Bridgman measurements, are limited to 10 GPa and 
it would be interesting to increase the pressures to at least 
15 GPa. We have also fitted the Murnaghan equation of 
state to rocksalt structure data below 15 GPa: 

" m(^ + lV 
y = Voexp ^ L (1) 

achieving good fit for Ai=41.01 GPa and ^ 2 =4.35. 

As seen from Fig. [5] the calculated relative volume 
change from NaCl to CsCl structure at the transition 
pressure P t =8.3 GPa is a 6% compression. We are not 
aware of any experimental measurements of this value, 
however, we can estimate it from the Bridgman data£ to 
be at least 3%. 



III. ELASTIC CONSTANTS 

The behaviour of elastic properties of the AgBr is spe- 
cially interesting in the vicinity of the phase transition. 
Fig.[3]shows pressure dependence of the bulk modulus as 
calculated from the equation of state, as well as from the 
following relationship with elastic constants: 

B = -V^±= 1 -{c ll +2c l2 ) (2) 

The difference in the bulk modulus B between NaCl and 
CsCl structures is quite small. On the same figure [3] 
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FIG. 3. Bulk modulus of AgBr calculated for NaCl structure 
from P(V) function (solid line) and dj data (+) from Fig. 2] 
The bulk modulus for CsCl structure is calculated from P(V) 
function (dashed line) and from Cjj data (x) as well. Exper- 
imental data for T=195 K (o) and T=300 K (□) are taken 
from Loje^. The dotted lines denote experimental pressure 
derivatives of B at zero pressure. 

we have included two dotted lines showing experimen- 
tal pressure derivatives of the bulk modulus at P=0 and 
T=195 K and T=300 K.i From these data we see that 
our calculations reproduce the experimental bulk modu- 
lus quite well and that the results from different methods 
are consistent within 10% . 

We have used the standard stress-strain relationship to 
derive all three elastic constants as a function of pressure. 
The method is based on constructing a set of linear equa- 
tions from stress-strain relationships for several deforma- 
tions of the unit cell. This set of equations represents a 
general form of the Hook's law and can be solved with 
respect to the elastic constants. Since in practice this 
set of equations is over-determined, to solve it we have 
used a singular value decomposition algorithm which au- 
tomatically provides a least squares solution of the set^S. 
For each pressure three negative and three positive defor- 
mations of each kind, namely, bulk, shear and tetragonal 
elongation, have been set up. 

TableUcompares the calculated and experimental elas- 
tic constants at ambient pressure and two temperatures: 
T=195 K and T=300 K. Present results show quite good 
agreement with the measurements, although the calcu- 



TABLE I. Experimental and calculated bulk modulus B and 
elastic constants Cij of the AgBr crystal, in GPa. Values in 
parenthesis are calculated from B = ^(cii + 2ciz). Experi- 
mental data for T=195 K and T=300 K are taken from Loje 4 . 
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FIG. 4. Elastic constants en, C12, C44 of AgBr in NaCl (solid 
lines and o, □, o) and CsCl (dashed lines and A, <, V) struc- 
tures respectively. The lines are cubic fits to the calculation 
data. The vertical solid arrow denotes predicted shear insta- 
bility pressure. 

lations are performed at T=0. The temperature data 
points allow to extrapolate the experimental values to 
T=0. Thus, we have estimated that from T=195 K to 
T=0 K the elastic constants might increase by about 5%. 
The mechanism of the phase transition from NaCl to 
CsCl structure involves variation of the C44 elastic con- 
stant. We have calculated C44 for both NaCl and CsCl 
structures and for several pressures in the interval from 
to 40 GPa . The results are shown at Fig. @] We can 
clearly see that the two lines representing the C44 elastic 
constants , intersect at pressure P=21 GPa, and more- 
over at value close to zero. The pressure P=21 GPa of 
the shear instability and P=27 GPa of the free energy 
equality are not exactly the same. However, we shall see 
in the next section that in this pressure interval another 
phase could be more stable than the cubic NaCl or CsCl 
structures. In previous sections the calculation were per- 
formed for NaCl and CsCl structures under constrains 




FIG. 5. Relation between cubic NaCl and monoclinic KOH 
(left) and between cubic CsCl and monoclinic KOH (right) 
unit cells. The KOH unit cell is drawn with dashed lines. 
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FIG. 6. Transition from NaCl to CsCl structure. Solid lines 
are guide to the eye. 



that both phases remain always cubic. Now, we release 
this restrictions and allow the unit cell to deform to a 
structure of PI symmetry being intermediate configura- 
tion between the NaCl and CsCl ones. 



IV. INTERMEDIATE PHASE 

When the intermediate configuration becomes mon- 
oclinic it resembles the KOH-type structure. The re- 
lationship between NaCl, KOH and CsCl structures is 
shown on Fig(5] Setting the monoclinic angle (3—90° or 
/3=125,27° one reestablishes the NaCl or CsCl structures, 
respectively. Hull and and Keeni^ pointed out that this 
monoclinic phase, may provide a proper transition path 
between NaCl and CsCl structures. A detailed study 
of the above phase transition mechanism has been per- 
formed and we refer the results below. We have prepared 
a model of the AgBr crystal in the KOH-type structure 
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FIG. 8. Free energy difference in the transition of AgBr 
from NaCl (squares) to CsCl structure (circles) through KOH 
structure (crosses). The polynomial fit to the CsCl data 
points is taken as the reference line (zero level). The lines 
are guide to the eye. 



with cell volume corresponding to the NaCl structure at 
P=0. The f3 angle of the monoclinic unit cell was set 
to 110°. Then the crystal geometry inducing atomic re- 
laxations was optimised with respect to free energy for 
various increasing pressures from to 45 GPa. At each 
pressure the minimisation procedure was reiterated sev- 
eral times since large changes in the unit cell size in- 
volved changes of the Brillouin zone grid. The results 
are plotted in Fig [6l One notices that at about 8 GPa, 
there is a clear tendency to depart from the pure NaCl 
structure to an intermediate monoclinic structure. There 
is also a clear jump of lattice parameters at the P=25- 
30 GPa range. Simultaneously, the a, and 7 angles do 
not change by more then 0.05° for the scanned pressure 
range. To verify the last findings we have carried out 
two series of calculations each starting from the same ei- 
ther low or high pressure configurations. The starting 
low and high pressure configurations were taken from 
optimised monoclinic configurations at P; c =25 GPa or 
P/ii=30 GPa, respectively. The results are shown in 
Fig. [7] It is clear from this figure that both monoclinic 
structures may show weak metastable behaviour. How- 
ever, the energy difference is small and it is possible that 
this effect will manifest itself only in low temperatures. 
Furthermore, the energy of the high-pressure configura- 
tion is always above the low-pressure one. Consequently, 
the high-pressure structure is probably never present in 
the real crystal, and we will omit this structure from now 



TABLE II. Calculated elastic constants (in GPa) of the mon- 
oclinic, KOH-type structure of the AgBr under pressure of 
P=20 GPa. 



FIG. 7. Cell parameters in NaCl structure to CsCl structure 
transition of AgBr. Solid lines are fourth degree fits to the 
data points. 
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TABLE III. Calculated structural parameters of the intermediate KOH-type structure of the AgBr crystal. 



P [GPa] 


a [A] 


6 [A] 


c[A] 


p n 


b/a 


c/a 





4.0803 


4.1064 


5.7601 


90.0 


1.006 


1.412 


5 


3.9521 


3.9819 


5.5982 


89.9 


1.008 


1.417 


10 


3.8226 


3.9411 


5.4886 


91.6 


1.031 


1.436 


15 


3.7368 


3.9091 


5.4100 


96.1 


1.046 


1.448 


20 


3.5445 


4.0156 


5.3589 


101.4 


1.133 


1.512 


25 


3.4033 


4.1011 


5.3403 


107.0 


1.205 


1.569 


30 


3.2993 


4.1976 


6.7106 


132.7 


1.272 


2.034 


35 


3.2316 


4.2612 


6.2892 


130.4 


1.319 


1.946 


40 


3.1725 


4.4554 


5.5465 


125.6 


1.404 


1.748 


45 


3.1489 


4.4375 


5.4855 


125.5 


1.409 


1.742 



on. Difference of the free energy of the three optimised 
candidate structures (NaCl, KOH, CsCl) are shown in 
Fig. [8j The polynomial fit to the CsCl free energy data 
is taken as the reference line ("0" level). At zero pres- 
sure the NaCl structure is the most stable one. When 
we increase the pressure above about 8 GPa, the cu- 
bic symmetry is broken, the monoclinic angle (3 starts 
to increase and the monoclinic phase is favoured. At 
the special point at P=21 GPa, at which the shear elas- 
tic constants C44 of NaCl and CsCl structure vanish on 
the Fig. |31 the monoclinic phase is still the most stable 
one. Only above P=35 GPa the CsCl structure becomes 
the most stable phase. However, for the finite temper- 
ature case this might change, since the NaCl and CsCl 
structures have different number of atoms in coordina- 
tion shells and hence different entropy terms in the Gibbs 
free energies. Since there are no low temperature mea- 
surements of pressure dependent behaviour of the AgBr 
we cannot compare the present results with experimental 
findings. 

We have also, for completeness, calculated some elastic 
constants of the intermediate KOH phase of the crystal. 
The monoclinic phase has 13 independent elastic con- 
stants. However, due to the small energy difference be- 
tween phases connected by the shear deformation, only 
6 of them can be calculated accurately. The remaining 7 
constants are too small and could be only estimated to be 
below 10 GPa. This is quite expected result. The 6 large 
elastic constants, calculated under pressure of P=20 GPa 
are presented in the Table [TTJ The values are consistent 
with the values of elastic constants of NaCl and CsCl 
structures depicted on the Fig. @J 

Summarising, our results support the hypothesis for- 
mulated by Hull and Keen,— that the phase transition 
from NaCl to CsCl phase, goes through the intermediate 
KOH-type, monoclinic structure. 

V. CONCLUSIONS 

The present DFT calculations show that at tempera- 
ture T=0 K the cubic NaCl and CsCl-type structures of 
AgBr are stable at pressures 0-8 GPa and above 35 GPa, 



respectively. Moreover, the calculation performed in con- 
strains of these two cubic structures permit to find the 
equation of state and the elastic constants in the whole 
interesting range of pressures from to 45 GPa. It is 
interesting to see that the shear elastic constants C44 for 
both phases NaCl and CsCl vanish close to 21 GPa. The 
derived equation of state and the elastic constants charac- 
terise properly the stable cubic phase, and additionally 
can provide an useful estimations of similar quantities 
within the intermediate monoclinic phase. The calcula- 
tion of the elastic constants of the monoclinic KOH phase 
showed that, due to the small energy difference between 
phases, the KOH phase is very soft to shear deformations 
and its elastic constants are consistent with the elastic 
constants of the other two phases. 

The study of the intermediate phase has been done 
without point symmetry elements constraints. This 
means that any deformation of the simulated supercell 
could have taken place. The only imposed limitation was 
the number of atoms in the supercell. Due to this ap- 
proach the supercell remained always monoclinic (or cu- 
bic for special values of monoclinic angle /3). The values 
of the free energy of the three candidate phases indicate 
that the monoclinic phase becomes the most stable one 
in the pressure interval from 8 to 35 GPa. This struc- 
ture has previously been proposed to be the transfor- 
mation path between NaCl and CsCl-type phases. And 
indeed, the free energy of the intermediate monoclinic 
phase seems to be lower then these of the cubic ones. 
However the difference in energy is quite small, which 
means that defects and impurities present in the crystal 
could show a large impact on the monoclinic phase. 



ACKNOWLEDGMENTS 

The use of facilities of the ACC "Cyfronet", Cracow, 
where the calculations have been done, are kindly ac- 
knowledged. This work was partially supported by the 
Polish State Committee of Scientific Research (KBN), 
grant No 2 P03B 004 14, and computational grant No. 
KBN/SGI_OPJGIN_2000 /IF J / 128/ 1998. 



6 



* Pawel.Jochym@ifj.edu.pl 

1 R.K. Singh, P. Khare, phys. stat. sol. (b) 103, 337 (1981) 

2 D.C. Gupta and R.K. Singh, Phys. Rev. B, 43 11185 
(1991). 

3 J.R. Chelikowsky, Phys. Rev. B 35, 1174 (1987) 

4 K.F. Loje and D.E. Schuele, J. Phys. Chem. Solids, 31 
2051 (1970). 

5 P.W. Bridgman, Proc. Am. Acad. Arts Sci., 74, 21 (1940) 

6 P.W. Bridgman, Proc. Am. Acad. Arts Sci., 76, 1 (1945) 

7 T.E. Slykhouse and H.G. Drickamer, J. Phys. Chem. Solids 
7, 207 (1958) 

8 W.A. Bassett and T. Takahashi, Am. Mineral. 50, 1576 
(1965) 

9 B.M. Riggleman and H.G. Drickamer, J. Chem. Phys. 38, 
2721 (1963) 

S.N. Vaidya and G.C. Kennedy, J. Phys. Chem. Solids, 32 
951 (1971) 

1 T.A. Fjeldly and R.C. Hanson, Phys. Rev. B 10 3569 
(1974) 

2 B.E. Mellander, Phys. Rev. B 26 5886 (1982) 

3 S. Ves, D. Glotzel, M. Cardona and H. Overhof, Phys. Rev. 
B 24 3073 (1981) 

4 G.S. Nunes, P.B. Allen, J.L. Martins, Phys. Rev. B 57, 
5098 (1998) 



15 S. Hull, D.A. Keen, Phys. Rev. B 59, 750 (1999) 

16 K. Kusaba, Y. Syono, T. Kikegawa, O. Shimomura, J. 
Phys. Chem. Solids 56, 751 (1995) 

17 R.N. Schock and J.C. Jamieson, J. Phys. Chem. Solids 30, 
1527 (1969) 

18 J.S. Lin, A. Qteish, M.C. Payne and V. Heine, Phys. Rev. 
B 47, 4174 (1993). 

19 J. Goniakowski, J.M. Holender, L.N. Kantorovich, M.J. 
Gillan and J.A. White, Phys. Rev. B 53, 957 (1996). 

20 M.C. Payene, M.P. Teter, D.C. Allan, T.A. Arias, J.D. 
Joannopoulos, Rev. of Mod. Phys. 64, 1045 (1992) 

21 CASTEP is a part of Cerius 2 system developed by Molec- 
ular Simulations. 

22 G. Kresse and J. Furthmiiller, VASP program, Vienna 
(1999); Phys. Rev. B 54, 11169 (1996); Comput. Mat. Sci- 
ence 6, 15 (1996). 

23 H.J. Monkhorst and J.D. Pack, Phys. Rev. B 13, 5188 
(1976). 

24 S. Froyen, M.L. Cohen, J. of Phys. C: Solid State Physics, 
19, 2623 (1986). 

25 S. Froyen, M.L. Cohen, Phys. Rev. 29, 3770 (1984). 

26 W. H. Press, B.P Flamery, S.A. Tenkolsky, and W.T. Vet- 
terling, Numerical Recipes (Cambridge Univ. Press, 1988), 
p.301 



